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Abstract 

In the context of a Black-Scholes economy and with a no-arbitrage argument, 
we derive arbitrarily accurate lower and upper bounds for the value of European 
options on a stock paying a discrete dividend. Setting the option price error 
below the smallest monetary unity, both bounds coincide, and we obtain the 
exact value of the option. 



1 Introduction 

In the seminal paper of Black and Scholes (1973), the problem of valuing a European 
option was solved in closed form. Among other things, their result assumes that the 
stochastic process associated to the underlying asset is a geometric Brownian motion, 
not allowing for the payment of discrete dividends. Yet the majority of stocks on 
which options trade do pay dividends. 

Merton (1973) was the first to relax the no-dividend assumption, allowing for a 
deterministic dividend yield. In this case, he showed that European options can be 
priced in the context of a Black-Scholes economy, with either a continuous dividend 
yield or a discrete dividend proportional to the stock price. However, when the 
dividend process is discrete and does not depend on the stock level, the simplicity of 
the Black-Scholes model breaks down. 

Let St denote the value of the underlying asset at time t, and let T be the maturity 
time of the option. When the risky asset pays a dividend D at time r < T, a jump 
of size D in the value process happens at that point in time. The stock price process 
is discontinuous at i = r and is no more a geometric Brownian motion in the time 
interval [0,T]. 

The standard approximation procedure for valuing European options written on 
such a risky asset, first informally suggested by Black (1975), considers a Black-Scholes 
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formula, where the initial price of the underlying stock So is replaced by its actual 
value less the present value {PV) of the dividends (Div), 



So 



S^ = So- PF(Div) 



This adjustment is made to evaluate the option at any point in time before r. After the 

payment of dividends, there is no need for further adjiistments. In this approximation, 
the input in the Black-Scholes formula is the value of the (continuous) stochastic 
process, 



where r is the risk-free rate. 

For t < T, the discontinuous stock price process St can thus be seen as the sum 
of two components {St = St + -De^''^^^*'). One riskless component, _De^'"('^^*\ 
corresponding to the known dividends during the life of the option, and a continuous 
risky component St ■ At any given time before r, the riskless component is the present 
value of the dividend discounted at the present at the risk-free rate. For any time 
after t until the time the option matures, the dividend will have been paid and the 
riskless component will no longer exist. We thus have St = S^ and, as pointed out 
by Roll (1977), the usual Black-Scholes formula is correct to evaluate the option only 
if St follows a geometric Brownian motion. In that case, we would use in the Black- 
Scholes formula Sq for the initial value, together with the volatility of the process S^, 
followed by the risky component of the underlying asset. 

If we assume that 5*^* follows a geometric Brownian motion, a simple application of 
Ito Lemma shows that the original stock price process Si does not follow a geometric 
Brownian motion in the time interval [0,t[. On the other hand, mider the Black- 
Scholes assumption that St follows a geometric Brownian motion in [0, r[, the risky 
component S^ follows a continuous process that is not a geometric Brownian motion 
in [0, t[. Therefore, the standard procedure described above must be seen as an 
approximation to the true value of such calls under the Black-Scholes assumption. As 
argued by Bos and Vandermark (2002), this assumption is typically underlying the 
intuition of traders, but the approximation is sometimes bad. In fact, as noticed in 
the early papers about option pricing (Cox and Ross, 1976; Merton, 1976a; Merton, 
1976b) , the correct specification of the stochastic process followed by the value of the 
underlying stock is of prime importance in option valuation. 

The deficiency of this standard procedure is reported in Beneder and Vorst (2001). 
Using Monte Carlo simulation methods, these authors calculate the values of call op- 
tions under the Black-Scholes assumption, and compare them with the values obtained 
with the approach just described. Reported errors are up to 9.4%. They also find 
that the standard procedure above usually undervalues the options. For these rea- 
sons, Beneder and Vorst (2001) propose a difi^erent approximation, trying to improve 
the standard procedure by adjusting the volatility of the underlying asset. This ap- 
proach consists in modifying the variance of the returns by a weighted average of an 
adjusted and an unadjusted variance, where the weighting depends on the time r of 
the dividend payment. Performing much better than the former approximation, this 
method still does not allow the control of the errors committed for the given param- 
eters of the economy. Analogously, Frishling (2002) warns on the mispricing risk due 
to the use of an incorrect underlying stochastic process. This discussion is followed 
by a series of recent papers suggesting different approximations that better match 
numerical results (Bos and Vandermark, 2002; Bos et al, 2003). More recently, Haug 
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et al (2003) discuss this problem. However, as these authors claim, "[i]n the case of 
European options, the above techniques are ad hoc, but the job gets done (in most 
cases) when the corrections are properly carried out" . 

The development of these approximations enhance two important aspects. First, 
they are not exact, and it is not possible to control the error with respect to the 
correct value of the option. Second, there are numerical procedures to estimate the 
value of these options, as for example, Monte-Carlo simulation methods. However, 
this method is time consuming and provides a convergence of statistical nature. 

The purpose of this paper is to derive a closed form for the exact value of Euro- 
pean options on a stock paying a discrete dividend, in the context of a Black-Scholes 
economy. We obtain an exact result and we need not to rely on ad hoc assumptions. 

This paper is organized as follows. In Section 2, an integral representation for the 
value of European options written on an asset paying a discrete dividend is obtained, 
and the convexity properties of the solutions of the Black-Scholes equation are de- 
rived. In section 3, we construct functional upper and lower bounds for the integral 
representation of the value of an option. These bounds follow from a convexity prop- 
erty of the solutions of the Black-Scholes equation. Theorem 13.41 is the main result 
of this paper and gives the algorithmic procedure to determine the price of European 
options on a stock paying a discrete dividend. In section 4, numerical examples are 
analyzed and we discuss the advantages of the proposed method. In section 5, we 
summarize the main conclusions of the paper. 

2 Valuation of European options on a stock paying 
a discrete dividend 

In this section, following a standard procedure to derive the Black-Scholes formula 
(Wilmott, 2000), we derive an integral representation for the value of a European 
option written on an asset paying a known discrete dividend. 

We consider a European call option with maturity time T and strike price K . This 
call option is written on an underlying asset with value St-, with stochastic differential 
equation, 

dSt = uStdt + aStdWt 

where jj, and a are the drift and volatility of the underlying asset. The quantity Wt is a 
continuous and normally distributed stochastic process with mean zero and variance 
t. Under these conditions, the underlying asset with value St follows a geometric 
Brownian motion. We also assume a risk-free asset with constant rate of return r. 

In the context of the Black-Scholes economy, the value V of an option is dependent 
of the time t and of the price of the underlying asset S. Under the absence of arbitrage 
opportunities (Wilmott, 2000; Bjork, 1998), it follows that V{S,t) obeys the Black- 
Scholes equation. 

The Black-Scholes equation is a quasi-linear parabolic partial differential equation, 
with S* > 0, and t > 0. To determine the solutions of the Black-Scholes equation, we 
introduce the new variables, 

f e = T-t 



log^+ (r-^) (T-t) 
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together with the new function Lp{x, 6) ~ e^'^'^ ^'^V{S, t). In the new coordinates H2.2f) . 
the Black-Scholes equation (|2.1|l becomes the diffusion equation, 
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(2.2) 



where x G R and 61 > 0. If e = 0, by ^T^, we have <p(a;,0) = V{S,T), and 
ip{x,T) = e^^V{S,0). Therefore, by (|2.2() . the forward solution in the time 9 of 
the diffusion equation relates with the backward solution in the time t of the Black- 
Scholes equation 1)2. The Black-Scholes problem for the price of a call option is to 
determine the option value at time t — whose value at maturity time T is, 

y (5, T) = max{0, 5 - if} (2.3) 

Therefore, due to the change of coordinates H2.2|l . the call option solution of the 
Black-Scholes equation (|2.1(l is equivalent to an initial value problem for the diffusion 
equation. 

Suppose now an initial data problem for the diffusion equation (|2.2(l . ip{x,d — 
0) — f{x). Under these conditions, the general solution of (|2.2|l is (FoUand, 1995), 
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2a^e 

and the solution of the Black-Scholes equation for a call option is, 

{x - y) 



dy (2.4) 



V {S, 0) = e-^^^ (x, T) = — == / y(e^ r)exp 

crV27rJ J_oo 



2ct2T 



dy (2.5) 



This integral can be easily calculated to obtain the usual Black-Scholes formula (Black 
and Scholes, 1973; Wilmott, 2000). 

For a dividend distribution at some time r g (0,T), the Black-Scholes formula is 
no longer true, since, during the life time of the option, the value of the underlying 
asset does not follow a geometric Brownian motion. However, if we take the time 
intervals, Ii = [0, t[ and I2 = [t,T], the value of the underlying asset follows a 
geometric Brownian motion in each interval Ii and I2, and, at time t — t, it has a 
jump equal to the dividend D. 

Before considering this case, we proceed with some properties of the solutions (|2.4|l 
and H2.5() of the diffusion and of the Black-Scholes equations. 

Definition 2.1. A real valued function f (x), with x Cz H, is convex if, for every 

Xi,X2 & R, 

<^(/(^l) + /(^2)) 

A simple property of convex functions is that, if the real- valued functions / and 
g are both convex, and g is increasing, then f{g{x)) is also convex. 

Proposition 2.2. Let f{x) the initial data function of a well-posed diffusion equation 
problem, and suppose that f{x) is non-negative and convex. Then, for fixed 9, the 
solution ip{x,6) of the diffusion equation is also convex. Moreover, if f{x) is an 
increasing function, then, for fixed 9, ip(x,9) is also increasing. 
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Proof. Suppose that the sohition H2.4|l of the diffusion equation (|2.2|l is well defined 
(FoUand, f995). By (|2.4ll . with z = y ~ x, we have, 



As, by hypothesis, f{x) is convex, then, for every z G R, 

{xi + z) + {X2 + z) 



f 



and, as f{x) is non-negative, 



= /(^ + ^^) <l[fiz + x^) + fiz + X2)] 




y ( ^^^^^ ,s] = — / f ( ^ +-^-4r-^ ] -TTTTl ] dz 



2 ; "V 2a20 

< ^[ip{x,,e) + ^ix2,e)] 

and so 1^9(2;, 0) is also convex. Assuming now that f{x) is increasing, we have that 
/(a^2) > fi^i): whenever X2 > xi. Then, for every z e R, we have, f{z + X2) > 
{z + xi), and, by (|2.4I) . the last assertion of the proposition follows. □ 

As (|2.3|) is a convex function in 5, Proposition 12.21 implies that the backward 
solution (|2.5() of the Black-Scholes equation H2.ll) is also a convex function. 

Suppose now that a dividend on the underlying asset is distributed at time t = t. 
We denote this dividend by D. According to the classical solution of the Black- 
Scholes equation (Wilmott, 2000), the price of the option just after the distribution 
of dividends at time t = t is, 

V {S+,t) = 5+ TV + - Ke-'''^'^-'''>N (d) (2.6) 

where, 

^_ \nS+-\nK+{r-^a^) {T - t) 



and denotes the value of the underlying asset just after the dividend distribution. 
The function Af(-) is the cumulative distribution function for the normal distribution 
with mean zero and unit variance. By Proposition 12.21 the function V[S^.t') is 
convex. Note that, the solution (|2.6|) is given by, y (S'+,t) = e^'''^'^^^^((>{x,T — t), 
and is directly calculated from (|2.5() and H2.3() . 

The approach taken here to value an option is equivalent (see, among others. 
Cox and Ross, 1976; Harrison and Krebs, 1979) to write this value at any point in 
time as the expected discounted payoff of the option at maturity T, under the so- 
called risk-neutral probability measure. Hence, knowing beforehand the amount to 
be distributed as dividend, the value of the option is not supposed to jump at r. In 
other words, the payment of known dividends D at a known point in time r does not 
affect the expectations at time r about the final payoff of the option at maturity T, 
and the value of the option is continuous at (Wilmott, 2000, pp. 129-131). Going 

^According to Wilmott, 2000, the jump condition on the asset price is known a priori, implying 
that there is no surprise in the fall of the stock price. Therefore, in order to avoid arbitrage oppor- 
tunities, the value of the option should not change across the dividend date. This is a no-arbitrage 
argument. 
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backward in time, the value of the underlying asset jumps from S+ to 5_ = 5+ + D, 
where 5_ is the value of the underlying asset just before the dividend distribution. 
As V{S+,t) — V{S-, t), by H2.6|l . the price of the option just before the distribution 
of dividends at time t — t is, 

where, 

- ln{S--D)-\nK+{r-^a'^){T-T) 



(2.8) 

In Fig. O we plot V{S+,t), V{S-,t) and V{S,T) as a function of 5. The 
functions V{S+,t), V{S-,t) and V{S,T) are convex. 




Figure 2.1: Option values V{S+,t), V{S^ , r) and V{S, T) as a function of the value S 
of the underlying asset. Parameter values are: fi ~ 0.01, a — 0.2, r = 0.03, K = 100, 
£) = 5, T = 1 and r = 0.5. 



To calculate the value of a call option as a function of the actual price (t = 0) 
of the underlying asset, we must introduce the change of coordinates H2.2|l into (|2.7|l 
and integrate as in H2.5|l . By H2.5|l and (|2.7(l . it follows that the time-zero value of a 
European option written on an asset paying dividend D at time t = r is given by, 



y(5,0)=e-''Xa;,r) 



V^[5_(y),T]exp 



dy (2.9) 



which has no simple representation in terms of tabulated functions. By Proposition 
lO V{S,{)) is also convex. 



3 Accurate bounds for ^(5,0) 

As it is difficult to determine a close form for the integral representation of the option's 
value H2.9() in terms of tabulated functions, to estimate the value V (S, 0), we use the 
convexity property of V{S^,t) and its asymptotic behavior as — ^ oo. 

Lemma 3.1. If K > 0, then, in the limit S- — s- oo, V{S-,t) is asymptotic to the line 
V = {S- - D) - Ke-'^'^-^'>, and V{S-,t) > {sJ- D) - ife-^'f^-^). 
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Proof. In the limit S- ^ oo, d oo, and N{d) 1. Hence, by (|T7|l . V{S-,t) is 
asymptotic to the hne Vi = {S- — D) — Xe"''*^^"'^-'. To prove the second part of the 
lemma, first note that, if Fi = {S- - D) - Ke-''^'^-^^ < 0, then S- < D + iCe-'^^^"^). 
As V{S^,t) is non-negative, if S- < D + Ke-'''^^-^\ then V{S^,t) > Vi. 

Suppose now that S- > D + Ke-'''-'^-''\ By hypothesis, we assume that there 
exists some S- = S such that, V{S, t) = (S - D) - Ke'''^^-^^ and ¥{8, r) > 0. By 
and we then have. 



,^ . N \d(S) + (tVT~^ - 1 - 



N [diS)] ~ 1 

As {S- — D) > Ke^^'-'^^'^\ from the equality above, we obtain, 



N [diS)] - 1 

Hence, 



-{S -D) = Ke-"-^' < (5_ - D) 



N 



d{S) + aVr^ < N [d{S)] 



which contradicts the fact that N{-) is a monotonically increasing function of the argu- 
ment. Therefore, the function V{S-,t) and the line Vi = {S--D)-Ke-''('^-''> do not 
intersect for finite S. As V{S-,t) is a continuous function of then V{S-,t) > Vi 
in all the range of S-, and the lemma is proved. □ 

To estimate the solution ()2.9|l of the Black-Scholes equation, we use ProDOsition l2.2l 
and Lemma lH.ll to construct integrable upper and lower bound functions of V{S-, r). 
This constructions proceeds as follows. 

Let us choose a fixed number S- = S* > D, and divide the interval [D, S*] into 
M > 1 smaller subintervals. The length of the subintervals is AS — {S* — D)/M, 
and their extreme points are denoted by, 

S^ = D + iAS, i^O,...,M 

As the function V{S- , t) is convex, in each subinterval, the function V{S- , t) is 
bounded from above by the chord that connects the points {Si, V{Si, r)) and (5^+1, V{Si+i, r)). 
We define the constants, 

a^ = -^^^ [V{S,, r) - r)] , z - 1, . . . , M 

where by (|2.7(l . V{So,t) — 0. Therefore, in each interval [Si^i,Si], the function 
V{S-,t) is bounded from above by the function = ai{S- — Si-i) + V{Si-i,T). 

Let us define the characteristic function of a set / as, Xii^) = 1, if a; S /, and 
Xi{x) = 0, otherwise. Then, the function V{S-,t) in the interval [D,S*] is ap- 
proached from above by the piecewise linear function, 

M 

y+(5_,r) = ^[a.(5--5.-i) + V(5,_i,r)]x[s._„S.](^-) (3.1) 

i=l 

To extend the bound of V{S^, r) to 5_ > S* , we introduce the function, 

V+{S-,r) = [{S--S*) + V{S*,r)]x[s',co){S-) (3.2) 
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By Proposition 12 . 21 and Lemma ITTl for 5*- > 5*, V^{S-, r) is the chord connect- 
ing the point {S* ,V{S* ,t)) to the point at infinity. Therefore, we have proved the 
following: 

Lemma 3.2. The function V{S-,t) has the upper bound, 

V{S^,r)<V+{S-,T) + V+{S-,T), if S->D 

where and are given by i)^. J|] and respectively, and the function {V^ + 

V^) is piecewise linear and non-negative. If S- < D, V{S-,t) — 0. 

The construction of a lower bound for H2.7|l follows the same line of reasoning. 

In each subinterval [Si-i,Si] C [D,S*], we can construct a linear function that 
bounds from below the function V{S^,t). Due to the convexity of V{S-,t), we 
construct the lower bound through the derivative of V{S- , t) at the middle point of 
each interval [Si^i,Si]. We then have, 

M 

V{- (5_,r) = (^»+^'^) - ^^h) + ^ (^^+^'^)] X[5._„5.l(^-) (3.3) 

i=l 

where, 

r (5_,r) ^ , = -+N(d + <jVT^) 

and d is given by (|2.8|) . 

To extend the lower bound of y (S*- , r) to S- > S*, we use Lemma lSTI to introduce 
the function, 

V^'iS^,T)=[iS^-D)-Ke-^<^^-^^]x[S',^)iS-) (3.4) 

By Lemma [3. II V"2^(5_,r) bounds from below V{S-,t). Therefore, we have: 

Lemma 3.3. The function V{S^,t) has the lower bound, 

V{S-,t) > V,-{S-,t) + V^-{S-,t), if S->D 

where V-[' and are given by J^H.!^) and {S.^^ , respectively, and the function (V-[' + 
V^) is piecewise linear and non-negative. If S- < D, V{S-,t) — 0. 

Finally, we can state our main result: 

Theorem 3.4. We consider the Black-Scholes equation together with the ter- 

minal condition Y2.3\) . We assume that K > Q and a dividend D > is payed at 
the time r with < t < T. Let S = S* > D be a fixed constant and let M > 1 be 
an integer. Then, the solution of the Black-Scholes equation with terminal condition 
1^2. 5^) has the following upper and lower bounds: 

M 

i=l 

+ SN{d*) + e^""^ [V{S*,t) - S*] N{d* - cr^/^) 
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and 

M 



1=1 

M 

+ e-'--^[y(5,+ .,r)-y'(5,+ i,r)5,+ 
+ SN{d*) - e-''^ (d + Ke-''^^-^^^ N {d* - a^) 



where, 



di 
d 
d* 



M 

log5-log^, + (r + ia^)r 

log(5 -D)-\ogK+{r + \!j^){T - r) 
log5-log5* + (r + ia2)^ 



y(5,T) = {S~D)N{d)-Ke-''^'^-'''^N{d-aVT~^) 



(jy/2Tr{T-T) cr^27r(T-r) (S* - £») 
^ [m,r)-T/(^,_i,r)] 



A, = N{d,^i)-N{d,) 

and N(-) is the cumulative distribution function for the normal distribution with mean 
zero and unit variance. 

Proof. By Lemmata H3.2II and 

yf (^„,T) + y2"(5_,r) <l/(^_,r) <y+(5_,r) + y2+(5-,^), if S->D 

Multiplying this inequality by the factors as in the integral l|2.9|l . and integrating, we 
obtain the estimates of the theorem. □ 
Note that, for S* > D fixed, limM^oo MiS,0) 7^ liniM^ooVg. ^,^(5,0). How- 
ever, if S* is large enough, both limits can be made arbitrarily close. Technically, this 
is due to the way the exponential term in l|2.5(l contributes to the integral. 



4 Calculating the price of a call option on a stock 
paying a discrete dividend 

Theorem ll-i.4l is the necessary tool to determine the price of a call option when the 
underlying asset pays a discrete known dividend before maturity time T. In fact, 
Theorem 13.41 asserts that we can always find upper and a lower bound functions for 
V{S, 0), and the bounding functions approach each other as we increase M and S* . 
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Figure 4.1: Bounds V^, ^(5, 0) and V^, j^j{S, 0) for V{S, 0), calculated from Theorem 
13.41 for several values of S* and M. In a) we have chosen S* — D+Ke^^^^^'^'' ~ 103.5. 
In b), S* = 2{D + Ke-''^'^-'')) = 207.0. Parameter values are: = 0.01, a = 0.2, 
r = 0.03, K = 100, D = 5, T = 1 and r = 0.5. 
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Table 1: Bounds V^.^]^j{S,0) and ¥^,^^{3,0) for ^(5,0), calculated from Theorem 
18.41 for several values of S* and M, and S = 110. The exact value ^(5*, 0) has 
been obtained by the numerical integration of p. 9(1 . The interval error e is given by 
(|2.9(l . Parameter values are the same as in Fig. 14.11 and we have chosen S* = D + 
j^^^r(T-r) ^ iQ2,.h, S* = 1.5{D + Ke-''^'^-^'>) = 155.3 and S* = 2{D + Ke-''^'^-^'>) = 
207.0. 



s 


s* 


M 




V{S,0) 






110 


103.5 


10 


11.24 


12.87 


15.41 


4.166 


110 


103.5 


50 


11.61 


12.87 


15.35 


3.739 


110 


103.5 


400 


11.63 


12.87 


15.35 


3.721 


110 


155.3 


10 


11.39 


12.87 


13.20 


1.807 


110 


155.3 


50 


12.79 


12.87 


12.88 


0.096 


110 


155.3 


200 


12.87 


12.87 


12.87 


0.006 


110 


155.3 


400 


12.87 


12.87 


12.87 


0.002 


110 


207.0 


10 


10.64 


12.87 


13.45 


2.813 


110 


207.0 


50 


12.72 


12.87 


12.89 


0.170 


110 


207.0 


200 


12.86 


12.87 


12.87 


0.011 


110 


207.0 


400 


12.87 


12.87 


12.87 


0.003 



To determine the price of the option, we first choose fixed values for the approxi- 
mation parameters S* and M. If V^l m('^^ ^) ^'^'^ ^s* hii^^ 0) differ too much within 
some fixed precision, we then increase S* and M. 

To analyze the convergence of the functional bounds and to the true price 
of a call option, we take, as an example, the parameters: fi — 0.01 (drift), a — 0.2 
(volatility), r = 0.03 (interest rate), K = 100 (strike price), D = 5 (dividend), T = 1 
(expiration time) and r = 0.5 (time of dividend paying). In Fig. 14.11 we show 
Vjt j^.j{S,0) and Vg, ]^j{S,0), for several values of S* and M, and calculated from 
Theorem 13.41 Increasing M and 5**, the upper and lower bounds Vgl ^^(5', 0) and 
V^t ]^j{S,0) approach each other, increasing the accuracy to which the functionals 
bounds approach the option price. To quantify this approximation to the value of the 
option, we define the interval error as, 

e = \V+^MiS,0)~VsZ^,,{S,0)\ (4.1) 

In Tabled] we compare the values of the upper and lower bounds Vgt j^jiS, 0) and 
Vg, j^{S,0), calculated from Theorem 13.41 with the exact value of V{S,0), obtained 
by the numerical integration of 1)2. 9|) . We show also the interval error e associated 
to both bounds. Assuming an interval error below the smallest unit of the monetary 
currency, for example, e < 10~^, we obtain the true value of the option. Therefore, 
for a choice of 5** and M such that e < 10~^, the difference between Vgl jv^(<S', 0) and 
V^t ]^j{S,0), is below the smallest unit of the monetary currency, and the rounded 
values of Vgl j(,f (5*, 0) and Vg, m C*^' ^) coincide. This rounded value is the option value 
within the chosen monetary accuracy 
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To analyze the global convergence behavior of V^t ^^(5, 0) and Vg. ^^(5*, 0), we 
chose a fixed value of S, and we change the approximation parameters S* and M. In 
Fig. 14.21 we show Vgt 0) ^^'^ ^s* mC"^' ^ function of 5*, for several values 
of AI. Increasing M, the upper and lower bounds of V{S, 0) become close in a region 
of the S* axis. A choice of 5"* in this region, gives better bounds to the value of the 
option, for lower values of M (Table ^E^nd Fig. I4.2|l . 

For all the examples we have analyzed, a good compromise to determine the value 
of the caU option is to choose S* = 2{D + JCe-'^C^-^)). Then, increasing iW, the 
interval error decreases. Due to the fast computational convergence of the expressions 
in Theorem 13.41 bounds with interval error below the smallest unit of the monetary 
currency are straightforwardly obtained. 




100 150 200 250 300 350 400 



Figure 4.2: Bounds V^, j^{S, 0) and Vg, j^.j[S, 0) as a function of S* , for 5 = 110 and 
several values of M. The parameter values are the same as in Fig. 14.11 and Tabled 



5 Concluding remarks 

We have obtained an upper and a lower bound for the exact value of a call option on a 
stock paying a known discrete dividend at a known future time. We have assumed the 
context of a Black-Scholes economy, where, away from the dividend time paying, the 
underlying asset price follows a geometric Brownian motion type stochastic process. 
The upper and lower bounds both approach the exact value of the option when two 
parameters are varied. In practical terms, one of these parameters {S*) can be fixed to 
the value, 5* = 2 (D + Ke~^'^^~'^^) , where K is the strike, D is the dividend, r is the 
time of paying the discrete dividend, and T is the length of the contract. Increasing 
the second parameter M, we obtain bounds for the option value with increasing 
accuracy. If this accuracy is below the smallest unit of the monetary currency, both 
bounds coincide, and we obtain the exact value of the option. 

The technique used to construct these bounds relies on the convexity properties 
of the option value at maturity, and on a property of the Black-Scholes and diffusion 
equations that preserves the convexity of propagated initial conditions. Under this 
framework, a similar methodology can be used to determine the value of a put option 
on a stock paying a known discrete dividend at a known future time. 
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Prom the numerical point of view, the technique developed here reduces to the 

sum of a few Black-Scholcs type terms, whereas mimcrical Monte Carlo methods rely 
on the poor convergence properties determined by the classical central limit theorem. 
In our numerical tests for the determination of the exact price of a call option, the 
computing time of our technique (using the Mathematica programming language) 
is several orders of magnitude faster than the computing time of finite diferences 
integration algorithms and of Monte Carlo methods. 
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